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Abstract 

From the integration of non-symmetrical hyperboles, a one-parameter 
generalization of the logarithmic function is obtained. Inverting this func- 
tion, one obtains the generalized exponential function. We show that func- 
tions characterizing complex systems can be conveniently written in terms 
of this generalization of the exponential function. The gamma function is 
then generalized and we generalize the factorial operation. Also a very re- 
liable rank distribution can be conveniently described by the generalized 
exponential function. Finally, we turn the attention to the generaliza- 
tion of one- and two-tail stretched exponential functions. One obtains, 
as particular cases, the generalized error function, the Zipf-Mandelbrot 
probability density function (pdf), the generalized gaussian and Laplace 
pdf. One can also obtain analytically their cumulative functions and mo- 
ments. 
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1 Introduction 



The convenience of generalizing the logarithmic function has attracted the at- 
tention of researchers since long ago }31 ! and particularly in the last years [9l 
[231 Uni 123 ISZl [S]. In Physics, several one-parameter generalizations of the 
logarithmic function have been proposed in different contexts such as: non- 
extensive statistical mechanics [SHI [23 [HI [HI [IH]j relativistic statistical me- 
chanics jm [5Dj and quantum group theory [T]. Also, more sophisticated such 
as two-parameter |22| and three-parameter [21j generalizations have been pro- 
posed, each one including previous situations as particular cases. Examples of 
the convenience of these generalizations have been seen in different fields, as 
for instance: psychophysics [33, neuroeconomics [3S1[I3], econophysics [5], 
complex networks [TSJ[3H], population dynamics [231 [5S] etc. 

Here, our main objective is to show that the generalized the stretched ex- 
ponential function, written as probability density function (pdf), is convenient 
to describe complex systems. In Sec. [21 we show that from the integration of 
non-symmetrical hyperboles, one obtains a one-parameter generalization of the 
logarithmic function, which we call g-logarithm. This generalization coincides to 
the one obtained in the context of non-extensive thermostatistics '3T . Inverting 
the q-logarithm, one obtains the generalized exponential (q-exponential) func- 
tion. Some properties of these generalized functions are presented. In Sec. [3l 
the gamma function is generalized and this allows us to generalize the factorial 
operation, which does not match the previous definition ^32i i4i. but as we show, 
it is self consistent. In Sec. [H the very reliable rank distribution obtained by 
Naumis and Cocho [28' is conviently described by the g-exponential function 
stressing finite size effects. In Sec. [5l we first show that the Zipf-Mandelbrot 
function, which is a fingerprint of complex systems, can be conveniently written 
in terms of the q-exponential. Raising the g-exponential argument to a given 
power, one obtains a function that generalizes the stretched exponential func- 
tion. Its generating differential equation is then presented. In Sec. [6l we con- 
sider the pdf 's for continuous variables. First we consider the one-tail stretched 
exponential generalization and obtain analytically its cumulative function and 
moments. One obtains the generalized error function as a particular case and 
the Zipf-Mandelbrot pdf as another. Next, we consider the two-tail generalized 
stretched exponential pdf and obtain analytically its cumulative function and 
moments. One has the generalized gaussian and the generalized laplace pdf as 
particular cases. The characteristic function is analitically calculated. Our final 
remarks are drawn in Sec. [71 

2 The g-generalized functions 

From the integration of non-symmetrical hyperboles, we obtain a one-parameter 
generalization of the logarithmic function, which coincides to the one obtained 
in the context of non-extensive thermostatistics [37] . Inverting this function, 
one obtains the generalized exponential function. Some properties of these gen- 



2 



eralized functions are presented. 



2.1 g-generalized logarithm function 

In the one-parameter generalization we address here, the q-logarithm function 
\nq{x) is defined as the value of the area underneath in the interval 

r dt - 1 

This is exactly the same function obtained from the non-extensive statistical 
mechanics context [36l |37] , which uses q = 1 — q. 

The usual natural logarithm (In x) is retrieved for q — and a linear function 
for q = 1. Scaling and deformation of the variable x are given by: lng(aa;'') = 
P\ni3q{a^^^x), so that for /3 = —1, one has hiq{a/x) = — \n^q{x/a) and for the 
particular case a = 1, lng(a;^^) — — ln_g(x). 



2.2 g-generalized exponential function 

The q-exponential function eq{x) is defined as the value, in such a way that 
the area underneath fq{t) = l/t^~^, in the interval t G [l,eg(x)], is x. This 
is the inverse of the g-logarithm function eq[lnq(a;)] — x — ln^[e^(a;)] and it is 
given by: 

f for < -1 

^'^""^ \ lim,~,^q(l + g'x)i/9' forgx>-l " ^ ' 

Notice that (1 + qxY^^ is real only if qx > —1. 

This is a non-negative function [eq{x) > 0, for all q and x] and 6^(0) = 1, 
independently of the q value. For q = 0, one retrieves the usual exponential 
function and for g = 1, a linear function. Notice that taking the surface 
underneath fl/t^~^ to be unitary, one generalizes the Euler's number: Cq = 
69(1) = i^ + qY^^- An interesting property is that [e^(a;)]° = Cq/aiax), meaning 
that the g-exponential argument scaling corresponds to a power of a different 
g-exponential function. For a — —1, one has that e_g(— x) — l/eq{x). 

The derivative of the g-exponential function with respect to x is: deq{kx) /dx — 
k [eq{kx)]^^^ ~ k eg/(i_g)[(l — q)kx], so that it is the solution of the following 
non-linear first order differential equation: dy{x) /dx — ky^~^{x), which is a par- 
ticular case of Bernoulli's differential equation: [d/dx+p{x)]y{x) — q{x)y^~'^{x), 
with p{x) = and q{x) = k. Notice that k has the dimension of y^ over the 
dimension of x. This means that it sets up a scale (inversion dimension of x) to 
the problem only if g = 0. One important application of the above differential 
equation concerns reaction kinetics [30j . 
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The successive derivatives of the g-exponential function with respect to x 
are written as: 



d"eq(x) 
dx" 



ri-l 



Wil -kq) 



.k=0 



i_„,- ^ (-g)"r(n - 1/g) 

r(-i/g) 



[eg (a;; 



1 1 — 



(3) 



where to get a closed form, one has used : nfe=o(l ~ ^'l) ~ (~9)"r(rt — 
l/g')/r(— l/g), which can be verified by induction. As far as we are aware, 
this is a new closed analytical form. 

The g-exponential function can be developed in Taylor's series around a; = 0: 



eq{x) 



E 



d"eq(a;) 



da;'' 



x=0 



^ ^ - {-qYT{n-l/q)x- 



where we have used that 6^(0) = 1 in Eq. [31 



3 Generalized Factorial 

The g-gamma function is defined as: 



Yi,{a) = / dt f'-^eqi-t) , (5) 
Jo 

with a > 0. For q ^ 0, one retrieves the standard gamma function 2 , since 
eo{—t) = e~* and 1/q — > oo. For q > 0, one has: Tq{a — ^ 0) = 1/a, Tq{l) = 

J^^^ dt eq{-t) = 1/(1 + q) etc. 

We point out that ^-generalized gamma function is closely related to the beta 
function Tq{a) = B(a, l/q + 1)/^''. From the properties of the beta function, 
one obtains the following recursion relations for the g-gamma function: 

which generalizes the factorial operation: 

l + g> + l) ' 

so that (a — > 0)g! — 1/a. This definition does not match the one previously 
proposed for integer values [32l|4j: n\q = Cq [J2k=i ^^^qi^)]- 

As stressed by Borges [TO], other possible generalizations of the exponential 
function are defined as one writes the generalization of the factorial in their 
Taylor expansion. Writing the Taylor expansion of the g-exponential function 
around a; = (Eq. [J]) as: 



eqix) 



oo 



n=0 ' * 
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we define the g-factorial as: 

, ^ -r(-l/g) n! ^ n(n-l)_^! 
(-g)"r(n - II q) 1 - q{n - 1) ' 

This definition of the generalization of the factorial is compatible with the 
generalization of the factorial obtained from the generalization of the gamma 
function if one write: 



^ I ^ n{n-l)q\ 

1 + g[n + sgn((7)] 



(9) 



where sgn(x) — x/\x\. 

For real and complex values of a, we propose the following definition for the 
generalized factorial function: 



^ I ^ a (a - 

1 + g[a + sgn((7)] 



(10) 



4 Beta-like Distribution 



Let us turn our attention to discrete random variables, the rank distribution 
in particular. We show that the rank distribution obtained by Naumis and 
Cocho [28J is conviently described by the g-exponential function and finite size 
effects are stressed. This rank distribution is very reliable since it has an under- 
lining microscopic model and has been validated by a wide range of experimental 
data [28]. 

To simultaneously fit the beginning, body and tail of experimental rank 
distributions of complex systems, Naumis and Cocho [28] consider N inde- 
pendent subsystems with a large number of internal states. The rank r of a 
system property dependent on the internal states of the subsystems decays as 
a two-free-parameter beta-like function: /(?■) — K[\ — r/{R + l)]^r°', where 
K = 1/Ef=i[l - r/{R+l)fr-", with R < N is the maximal value of r, K is 
the normalization factor and the two free parameters are a and (3. As noticed 
by the autors, if i? > 1, then K « 1/B(1 -a,l + 13). 

Finite size effects are described by factor [l—r/{R-\-\)]^ « l—[3r/R « 6^''/'"°, 
for R^ r with ro = R/(3. In this way, we see that the rank distribution of f{r) 
is in fact a generalization of the standard technique of multiplying the power-law 
(r~") by an exponential cutoff: /(r) cx r~°'e~^^^° . 

If one writes: [1 - r/(i? + 1)]^* = [1 - (l//3)/3r/(i?-f 1)]'^ = e^/^[_/3r/(i?-|- 1)] 
then 

K 



f( \ ^ 



R+l 



R+l 



(11) 



meaning that the g-exponential function can properly take into account finite 
size effects. The exponential cutoff is retrieved where R and b grow, but R must 
grow faster than /3, to have R oo and /3 — > oo, but (3/R-^ Q. 
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5 Special Functions and Processes 

In what follows, we first show that the Zipf-Mandelbrot function, which is a 
fingerprint of complex systems, can be written in terms of the g-exponential. 
Next, raising the g-exponential argument to a given power, one obtains a func- 
tion that generalizes the stretched exponential function. Finally, we obtain the 
process (differential equation) of which the generalized stretched exponential is 
the solution. 

5.1 Zipf-Mandelbrot Function 

The envelope of a rank distribution of complex systems can sometimes be very 
well described by the Zipf-Mandelbrot function TT, which may be written in 
terms of ^-exponential function: 



with q = 1/7, a = c/7 and A = d/c^. 

We remark that Eq.[T2]has another interesting application in time-dependent 
luminescence spectroscopy. In this case the relaxation processes are known as 
Becquerel decay function ^ . 

5.2 Generalized Stretched Exponential Function 

In some instances, to model some characteristic a; of a complex system it is nec- 
essary a deformation in the argument of Eq. 1121 which leads to the generalized 
stretched exponential function: 



The usual stretched exponential function, also known as the Kohlrausch 
function [141 [8], is obtained from Eq.[T3]in the limit q ^ 0. Although the (usual) 
stretched exponential function has been used to describe relaxation processes 
in time-dependent luminescence spectroscopy j8], a generalization of the form 
g/3/a g-ei/^(2;/Q)/a gggjjjg ]-,g jnore Convenient to fit experimetal data [6 . In 
this case we stress that the (j-exponential figures as the argument of the usual 
exponential function. 

The stretched exponential function can also be obtained as the solution of the 
following differential equation: dy{x)/dx = —x^/^~^y{x)/{aj3), which can be 
written in terms of the logarithmic function as: d\ny{x) / dx = —x^/^~^/(al3). 
If one replaces the logarithmic function in the differential equation obtained 
above by the g-logarithm, one obtains: 






d hi_qy{x) 
dx 



X 



a/3 



(14) 
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or equivalently 

which is a particular case of Bernoulh's equation: dy{x)/dx + p{x)y{x) = 
q(x)y^^^(x), with p{x) = and q{x) = x^^^~^ /{afS) and which solution is 
precisely the generalized stretched exponential function of Eq. [131 



6 Probability Functions 

Considering the factor A > and using the non-negativeness of Eq. 1131 we write 
the generalization of the stretched exponential pdf and study its properties. We 
consider one- and two-tail distributions and obtain some known pdf 's (general- 
ized gaussian) and new ones (generalized error function and generalized Laplace 
pdf) as particular cases. 



6.1 One-Tail PDF 



If the considered independent variable x is constrained to non-negative (or even- 
tually non-positive) values, then one uses the one-tail pdf. We consider the 
generalization of the stretched exponential pdf and analytically obtain its cu- 
mulative function and moments. From this pdf, one obtains the generalized 
error function as a particular case and the Zipf-Mandelbrot pdf as another. 
The normalization factor A of Eq. [TSlis: 



1 



dt e_ 



B 



q 



(16) 



where B(a, b) = T{a)T{b)/r{a+b) is the beta function [2 and r(a;) = dt V-'^e 
X > is the gamma function [2j . The integral of Eq. [1^] docs not diverge only 
ifO<q<l,0</?<l/g, and a > and one has the generalized stretched 
exponential pdf : 



d(1) 
q,a.,(3 



(x) = 



[q/af 



/3B(/3,l/q-/3) 
For our purposes, it is more convenient to write a 



a(3E{(3,l/q-f3) 



(a/3)'' > 0: 



(17) 



(18) 



which is depicted in Figs. [T] and [H 

As g ^ 0, B(/3, l/q — (3) « ql^T{f3), since 1/q » f3, one retrieves the stretched 
exponential function Po^^j .(a;) = /3''e-^(^/'')'^''/[a/3r(/3)] or pI^^^Jx) = e'"''^^ /[af^ I3T{(3)] 



The cumulative function of Eg. [TSlis: 



dt P 



(1) 



q,a,l3 



it) 



A 
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qx^/P 
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where 

2Fi{a,b;c;x) = }^—- — (19) 

n— 

is the hypergeometric function 2 and 

{a + n-iy. r{a + n) 



(a)n = a(a + l)(a + 2) • • • (a + n - 1) 

is the Pochhammer symbol. 
The moments of x are: 



{a-iy. Via) 



)q,a,l3 



my 



B[(n+l)/3,l/g-(n+l)/3] 

B(/3,l/q-/3) ' ^ ^ 

where one sees they are finite only if g < + l)/3]. 

If g < 1/(3/?), the mean value and variance are finite and, respectively, given 

by: 

^"^-^^ - UJ B(/3,l/g-/3) 
^-.a./3 ^ B(3/?,l/g-3/3)B(/3,l/g-/?) ^ 
3^(2/3,1/9-2/3) 

Notice that the ration cr? ^ p/ix}^ a p depends only on /3 and q, but not on a. 

Particular values of (3 leads to a g-generalization of the error function and 
to the Zipf-Mandelbrot pdf. 

6.1.1 Generalized Error Function. 

To generalize the error function, consider (5=1/2 and a = (3^ (or a = 1) in 
Eg . [T9l and one has: 

- ^i.'.WW - B(l/2. l)?- 1/2) /'"'^-iH')- (^S) 

As(7^0,B(l/2, 1/g— 1/2) = \/qTT and one retrieves the standard error function: 
erfo(a;) = erf (a;) (2/0F) /q" dt e"*'. 

6.1.2 Zipf-Mandelbrot PDF. 

For /3 = 1 in Eq. [181 one obtains the Zipf-Mandelbrot 's pdf: 

= (— ) ^ (1 + 5~a;/«)i/9 ■ 

which mean value and variance are: 



o'ifv.i a-1 



(26) 
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which is finite for < 5 < 1/3 and, from Eq. (THl one obtains its cumulative 
function: 



1 



where the upper-tail distribution its simply given by: 



Fq,a,lix) 



dtP^:Ut) = [l + q- — 



(27) 



(28) 



which is a more suitable for fitting the model to real data than the pdf (Eq. [24] 
itself. 



6.2 Two Tail PDF 

If the domain of the considered independent variable is not bounded, it is inter- 
esting to consider its absolute value |a;| in Eq . [TSl and one has a symmetric pdf 
about the line x = 0. Notice that in this case, the normalization factor must 
be halved since the domain has been doubled in a symmetrical way, the pdf is 
then: 



m 



P 



•B{/3,l/q-(3) ' 
which is a convenient function to use in wavelets 



-f3 



1//3' 



(29) 



Its cumulative function 



is: 



1 



sgn(x)4^2_^(|x|) 



(30) 



where F^^^ ^(x) is given by Eq. ]T9\ 

On one hand, due to its symmetry around x ^ 0, the odd moments of this 
pdf vanish (a;^""'"^) = 0, with n — 0,1,2,.... On the other hand, the even 
moments (a;^") are finite only if g < l/[(2n -I- !)/?]: 



-I 2n 



B[(l + 2n)/3,l/g- {l + 2n)/3] 
B{(3,l/q-p) 



(31) 



In the following we retrieve the generalized gaussian as a particular case 
of Eq. [29l Also, we propose to consider another particular case of Eq. [29] to 
generalize the laplacian pdf. The characteristic function of both particular cases 
are analytically calculated. 



6.2.1 Generalized Gaussian. 

An interesting particular case of the two-tail generalized stretched exponential 
function is when (3 = 1/2, which leads to the g-guassian [55] : 



Gq^aix) 



Vqni/q) e_,-[-(x/a)V2] 



r(l/g-l/2) 



(32) 
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Due to symmetry, all odd moments vanish. Even moments are given by Eq. 1311 
, ^2a^V r(l/2 + n)r(l/g-l/2-n) 

j V^r(i/,~-i/2) ■ 



and the variance is finite only for q < 2/3: 



2a 



2 



4,a = {x')^,a = 2^ . (34) 

Using that for a > 6, T{a+b)/r{a) = a^ when 9 ^ in Eq.[32l r(l/q)/r(l/g- 
1/2) = and one has a guassian, with variance : Go.ai^) = e"'^/"^' /^/v2t 

as a particular case. Another particular case is when q = 1 and one retrieves 
the Lorentzian (Cauchy pdf) Gi^a{x) = l/{7ra%/2[l + {x/a)'^/2]}. 

The characteristic function of Eq. [32] has an analytical closed form: 



Pg,a,i/2ik) = (e''^^)=/ dxQ,,(x)e^'=^ 

^ — OC 

[ 2 j r(i/g-i/2) 

where K^{z) is the AT-modified Bessel function [5]: 



(35) 



_ 2-r(. + i/2) dtcost 

For the Gaussian one has also a Gaussian po,a{k) = e"'^^'^^/^ but for the 
Lorentzian, one has the Laplace function pi^a — e^"'''!. 

6.2.2 Generalized Laplace PDF. 

For /? = 1 and g = 0, Eq.Hileads to the Laplace pdf P^'^l^ix) = e"l^l/''/(2a). 
For arbitrary q, one has the generalized Laplace pdf: 

^ PilM - ^ (^^) , (37) 

and its cumulative function is: 

F^lii^)-lhaM), (38) 

where Pq,a,i{x) is given by Eq. [Ml 

The odd moments of Eq. [37] vanish and the even ones are finite if g < 
1/(1 + 2n): 

{x^n4,a,i = ^^ {^) B(l + 2n,l/q-2n-l) . (39) 
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Its characteristic function has an analytical closed form: 

Pq,a,i{k) — I dx Lq^aix)e^''^ — / dx cos{kx) e-q{~ X / o) 



„2kx ^Ir^ 

-oo a 

1 - q f°° , cos{kx) 
dx 



(1 + qx/ay/'} 



^ty«^ Sin(to + 7r/(29)) , ^ / 13 1 , -<o 



. , ^ . , ^ , 1F2 1:1 , : — (40) 

sin(7r/^)r(l/^- 1) 1 ^ \^ ' 2g'2 2q i J ^ ^ 

where to = ka/q and iF2(a; 61, 62; z) = I]^=o{[(«)"]/[(^i)"(^2)«]}(a;")/n! is the 
hypergeometric function [2], with (a)„ being the Pochhammer symbol. 



7 Conclusion 

We have shown that the ^-generalization of the exponential is suitable to de- 
scribe complex systems since it generalizes functions, which govern these sys- 
tems. As special cases we have obtained the g-generalized stretched exponential 
function, which has the generalized error function, the generalized Laplace pdf 
and the already known generalized gaussians as special cases. Further, we have 
used the g-exponential to write the very reliable rank distribution obtained by 
Naumis and Cocho. An interesting application of the g-exponential is in the 
generalization of the gamma function, which allows us to obtain a new gener- 
alization for the factorial operation. Since these distributions are the solution 
of differential equations that describe the complex systems, the q generalization 
brings many different systems to be described by the same underlying process. 
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Figure 1: Behavior of Eq. [TH]as a function of: (a) x and q, with /? = 1/2 and 
a — 1, (b) X and /3, with q = 1/5 and a = 1, (c) q and /3, with x — 1 and a = 1. 




0.D 0,S l.ti 1.5 2.0 2,6 ^5 4.0 



(a) 

1,0 J , 




0,0 0.5 1,0 1.5 2,0 2.5 3.0 3,5 4,0 

sr 

(b) 

Figure 2: Behavior of Eq. [18] as a function of: (a) x with (3=1/2 and a = 1 
and for several q values (Projection of Fig. [1^), (b) x with q = 1/5 and a = 1 
and for several /? values (Projection of Fig. [Us). 
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